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ABSTRACT 

We measure equivalent widths (EW) - focussing on two unique features (Nal 
and Ti02) of low-mass stars (< O.3M ) - for luminous red galaxy spectra from 
the the Sloan Digital Sky Survey (SDSS) and X-Shooter Lens Survey (XLENS) 
in order to study the low-mass end of the initial mass function (IMF). We com- 
pare these EWs to those derived from simple stellar population models com- 
puted with different IMFs, ages, [a/Fe], and elemental abundances. We find that 
models are able to simultaneously reproduce the observed NaD A5895 and Na I 
A8190 features for lower-mass (~ cr*) early-type galaxies (ETGs) but deviate 
increasingly for more massive ETGs, due do strongly mismatching NaD EWs. 
The Ti02 A6230 and the Na I A8190 features together appear to be a powerful 
IMF diagnostic, with age and metallicity effects orthogonal to the effect of IMF. 
We find that both features correlate strongly with galaxy velocity dispersion. 
The XLENS ETG (SDSSJ0912+0029) and an SDSS ETG (SDSSJ0041-0914) 
appear to require both an extreme dwarf-rich IMF and a high sodium enhance- 
ment ([Na/Fe] = +0.4). In addition, lensing constraints on the total mass of 
the XLENS system within its Einstein radius limit a bottom-heavy IMF with 
a power-law slope to x < 3.0 at the 90% C.L. We conclude that Nal and TiO 
features, in comparison with state-of-the-art SSP models, suggest a mildly steep- 
ening IMF from Salpeter [dnjdm oc m~ x with x = 2.35) to x ~ 3.0 for ETGs in 
the range a = 200 — 335 km s -1 . 

Subject headings: dark matter — galaxies: elliptical and lenticular, cD — grav- 
itational lensing: strong — galaxies: kinematics and dynamics — galaxies: evo- 
lution — galaxies: structure 

1. Introduction 

When constraining the star formation, metallicity and gas/dust content of galaxies, 
the initial mass function (IMF) is often assumed to be universal and equal to that of the 
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solar neighborhood (Kroupa 2001; Chabrier 2003; Bastian, Covey & Meyer 2010). However, 
evidence has recently emerged that the IMF might evolve (Dave 2008; van Dokkum 2008) or 
depend on the stellar mass of the system (e.g. Worthey 1992; Trager et al. 2000b; Graves et 
al. 2009; Treu et al. 2010; Auger et al. 2010b; Napolitano 2010; van Dokkum & Conroy 2010). 
van Dokkum & Conroy (2010; hereafter vDCIO) suggested that low-mass stars (< 0.3 M ) 
could be more prevalent in massive early-type galaxies. The increase in the mass-to-light ratio 
(M/L) of galaxies with galaxy mass may thus be partly due to a changing IMF rather than 
an increasing dark matter fraction, consistent with previous suggestions (Treu et al. 2010, 
Auger et al. 2011, Barnabe et al. 2011, Dutton et al. 2012, Cappellari et al. 2012). vDCIO 
showed that some spectral features, such as the Na I AA8183, 8195 doublet (called NaI0.82 
by CvD12), depend strongly on surface gravity at fixed effective temperature, betraying the 
presence of faint M dwarfs in integrated light spectra. If correct, the low-mass end of the IMF 
can be inferred directly from red/near-IR spectra of old populations. Hence, the strength 
of the Na I doublet versus another sodium feature, such as the NaD doublet (called Na0.59 
by CvD12), should provide a powerful means for separating the IMF from other effects. 
Specifically for the purpose of determining the low-mass IMF down to ~ 0.1M o for metal-rich 
stellar populations with ages of 3-13.5 Gyr, Conroy & van Dokkum (2012; hereafter CvD12) 
presented new population synthesis models. The NaD feature responds more strongly to 
Na-enhancement than IMF in the CvD12 models, while the Na I doublet is strong in stars 
with mass < 0.3 M and weak or absent in all other types of stars. Unfortunately, NaI0.82 
is also sensitive to age and metallicity, and NaD is influenced by any interstellar medium. It 
is therefore necessary to test these models over a range of age and metallicity indicators, as 
well gainst other lines caused by low-mass stars. 

In this letter, we focus on the Nal feature as indicator of low-mass stars. We use NaD 
as indicator of a change in sodium abundance and H/3 and [MgFe] as indicators of age and 
metallicity, respectively. This allows us to assess model degeneracies and deficiencies. We 
propose the use of the TiO feature at A6230 as an indicator of the presence of low-mass 
stars. We find that both of these features (Nal and TiO) correlate with galaxy velocity 
dispersion, implying a steepening of the IMF slope in ETGs with a > a*. We assume 
H = 70kms _1 Mpc -1 , Q m = 0.3 and f2 A = 0.7 throughout this letter. 
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Fig. 1. — Galaxy (continuous lines) and model (dashed lines) spectra in the regions of the 
NaD (top) and Nal (bottom) features. The observed NaD EWs do not match the models for 
the most massive ETGs (> 300 km s -1 ). Na I absorption is stronger in the XLENS system 
and SDSSJ0041-0914 and appears in both cases to require an IMF steeper than Salpeter, 
while the stacked SDSS spectrum shows a weaker Na I feature that matches a model with a 
Salpeter IMF. The bottom panels show the noise spectrum of the XLENS system. 
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2. The data 

As part of the XLENSQ project, we obtained a UVB-VIS X-shooter spectrum of the 
massive and luminous early-type SLACS (Sloan Lens ACS Survey, Bolton et al. 2006) lens 
galaxy SDSS J0912+0029 at z = 0.1642, with high enough signal-to-noise to perform stellar 
population analyses. The lens galaxy shows a surprising ly deep NaI0.82 feature (FigQ, 
making it an extremely interesting target for studying the low-mass end of the IMF in 
ETGs. We measure the luminosity-weighted velocity dispersion of the lens galaxy from the 
reduced flux-calibrated ID UVB-VIS spectrum using the Penalized Pixel Fitting (pPXF) 
code of Cappellari & Emsellem (2004). We obtain (<7*)(< R cS ) = 325 ± 10 ± Ukms' 1 , in 
agreement with the previously published value (a ~ 313 ± 12kms _1 ; Bolton et al. 2006). 
We also used the spectra of ~250 galaxies with similar morphology and colors (all LRGs) 
from the Sloan Digital Sky Survey DR8 (SDSS; Aihara et al. 2011), in five velocity-dispersion 
bins spread over 200-335 km s -1 (~50 galaxies per bin). We examine one system, SDSSJ0041- 
0914, separately because it has a NaI0.82 feature comparably deep to the XLENS system. 



3. Stellar Population Synthesis Modeling 

We use the synthetic spectra of CvD12 to analyze the stellar populations of these galax- 
ies. The models make use of two separate empirical libraries, the MILES library covering 
3500-7400 A (Sanchez-Blazquez et al. 2006) and the IRTF library of cool stars covering 
8100-24000 A (Cushing et al. 2005; Rayner et al. 2009). They also incorporate synthetic 
spectra with the purpose of investigating changes in the overall metallicity or changes in 
the abundances of individual elements and to cover the gap in wavelength between the two 
empirical libraries. We refer to CvD12 for details. The abundance variations of single el- 
ements are implemented at fixed [Fe/H], which implies that the total metallicity Z varies 
from model to model. We measure line-strength indices in the range 4000-8400 A, including 
the standard Lick indices H/3, Mg6, Fe5270, Fe5335, NaD and a TiO index (Ti02) using the 
definitions of Trager et al. (1998), and the commonly-used [MgFe] combination^} 

We define a modified index around the Na I doublet 8183, 8195 A, which seems to be 
strongly dependent on the low- mass end of the IMF (Table 1). This index is slightly different 
from that used by vDCIO and CvD12, having a wider central index bandpass and slightly 
wider pseudo-continua. Our definition is more stable against velocity dispersion variations 
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and more suitable for massive ETGs. We convolve all the galaxy and model spectra to an 
effective velocity dispersion of a = 335 km s -1 (the upper limit in our sample), to correct 
for kinematic broadening, before measuring indices. Indices in both the observed and syn- 
thetic spectra are measured with the same definitions and method (SPINDEX2; Trager et 
al. 2008). We do not place our indices on the zero-point system of the Lick indices and quote 
them as equivalent widths (EWs) in units of A, except for Ti02, which is given in magnitudes. 



H/3 is primarily an age indicator, while a combination of Mg6, Fe5270, and Fe5335 
yields information on the mean metallicity [Z/H] of the population (Worthey 1994) while 
minimizing the effects of abundance ratio variations (e.g., Gonzalez 1993; Trager et al. 2000a). 
These indices (Panel (a), Fig. [2]) show a good agreement between the models and the galaxies 
EWs for old stellar populations, with an age of 13.5 ± 3 Gyr for a > 300 km s -1 (black 
points) and younger ages for lower mass ETGs. The statistical error is deduced directly from 
variations in H/3[^J The most massive ETGs have values of [a/Fe] between solar and super- 
solar (~0.2), in good agreement with the prediction that massive galaxies have significantly 
super-solar abundance ratios because of rapid, high-efficiency star formation (Trager et al. 
2000b; Thomas et al. 2005, Spolaor et al. 2009, 2010). Given the uncertainties in the line- 
strengths of the two individual galaxies (SDSSJ0912+0029 and SDSSJ0041-0914), we are 
unable to determine their ages and metallicities precisely, but their line strengths are similar 
to the mean of the highest-mass SDSS sample, with a deviation from the average EW smaller 
than la in both age and metallicity. The Nal and NaD indices can in principle be used to 



3 For stellar populations with ages > 10 Gyr, an uncertainty of 0.1 A in H/3 corresponds to 1 Gyr 
uncertainty in the age (cf. Worthey 1994). 



4. 



Results and discussion 



Table 1. Definition of the index around the Na I doublet 8183, 8195A 



Index 



Central band 



Pseudocontinua 



(A) 



(A) 



Nal 8168.500 - 8234.125 8150.000 - 8168.400 

8235.250 - 8250.000 
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Fig. 2. — Index-index plots of the main absorption features. Lines and crosses are different 
SSP models from CvD12 with increasing IMF (Chabrier, Salpeter with a slope of x = 2.35, 
a bottom-heavy IMF with slope of x = 3.0 and an extremely dwarf-rich IMF with a slope of 
x = 3.5). Points colored according to their velocity dispersions are individual SDSS galaxies, 
with index errors similar to SDSS J0041-0914. In the plots showing sodium, the XLENS 
system SDSS J09 12+0029 requires a very steep IMF, violating lensing contraints on its total 
mass (see the text for further details). Panel (a): H/3 as a function of [MgFe]. The most 
massive ETGs (> 300 km s -1 ) best match an old stellar population (13.5 Gyr) with super- 
solar total metallicity. Lower-mass systems are younger. Panel (b): NaD as a function of 
Nal. Only low-mass (< 250 km s -1 ) systems match the models. More massive ETGs require 
a higher [Na/Fe] and the XLENS system and SDSS J0041-0914 also require a very steep IMF 
slope. Panel (c): Ti02 as a function of [MgFe]. The most massive ETGs require an IMF 
slope slightly steeper than Salpeter. A Chabrier-type IMF systematically underestimates the 
SDSS Ti02 EWs. Panel (d): Ti02 as a function of Nal. The ETGs match with the models 
using a Salpeter or slightly steeper IMF, but the XLENS system and SDSSJ0041-0914 still 
do not match the SSP models well. 
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constrain the IMF slope (CvD12), and this relation is shown in Figure 2b. Although the data 
match the models for low-dispersion systems (<250kms _1 ), the models with solar [Na/Fe] 
abundance do not match the NaD strengths and only models with [Na/Fe] = +0.3 — +0.4 
dex match the NaD indices for higher-mass ETGs. We suggest two possible explanations for 
this behavior: 

(i) NaD is highly contaminated by the interstellar medium (ISM) for higher mass ETGs; 
for example, dust lanes provide additional absorption in this resonance line (Sparks et al. 
1997). Interstellar absorption within a galaxy may alter the stellar absorption profile and 
therefore the calculated EW, leading to an incorrect inference of the underlying stellar pop- 
ulation. 

(ii) Very massive ETGs have higher [Na/Fe] abundances (> 0.3 dex) and slightly 
bottom-heavy IMFs which correlate with their stellar velocity dispersions. Therefore, if 
we explain the strengths of these features in giant ETGs using abundance ratios, we require 
an average iron abundance in excess of solar ([Fe/H] ~ 0.2), a IMF with x = 3.0 and a 
high sodium abundance ([Na/H] > 0.3). In the a— enhanced bulge of the Galaxy, Fulbright 
et al. (2006) find an averaged [Na/Fe] =0.2 dex, and that [Na/Fe] < 0.3 dex in all stars. 
However, Lecureur et al. (2007) find that [Na/Fe] ratios increase sharply with metallicity. 
They obtain values of [Na/Fe] ~ 0.5 for [Fe/H] = and even higher for [Fe/H] > 0, but 
with a scatter of 0.29 dex resulting in a range of [Na/Fe] from —0.1 to almost 1.0. It is 
therefore possible for massive ellipticals to have high [Na/Fe]. In both cases, the models 
seem consistent with a Salpeter IMF at the low-dispersion end and a slightly bottom-heavy 
IMF for the high-dispersion end, if these effects are accounted for, but the models predict 
a steeper IMF slope of x ~ 3.0 - 3.5 for both the XLENS galaxy SDSSJ0912+0029 and 
SDSSJ0041-0914. 

We note that a TiO feature at 8 199 A could partly contaminate Nal, although this feature 
should not vary strongly (CvD12). To test possible contamination, we use a model with 
[Ti/Fe]=±0.3 and calculate the Nal EW for a Chabrier IMF. We find that Ti enhancement 
only affect the Nal index by 1%. 

Overall we conclude that the NaD EWs and its trend with stellar mass remain unex- 
plained for systems with a > 250 km s -1 . We find that SSP models predict that TiO features 
also depend strongly on the slope of the low-mass end of the IMF, such as Ti02, shown in 
Figures 2c and 2d. This indicator gives more support to the conclusion that the sodium 
strengths of the XLENS ETG, SDSSJ0912+0029, still remain somewhat difficult to explain 
by current stellar population models, although most SDSS systems can be matched in Nal 
for most ETGs (if not in NaD). Together the Ti02 and Nal indices both imply a bottom- 
heavy IMF, steepening from Salpeter to possibly x ~ 3 for the most massive SDSS ETGs. 
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As in Treu et al. (2010), a bottom-light IMF such as Chabrier IMF is inappropriate for the 
most massive ETGs. 



4.1. Limits on the IMF from Strong Lensing 

A strong case against an extreme bottom-heavy IMF can be made using the system 
with the strongest Nal EW (Fig. g), the XLENS galaxy SDSSJ0912+0029. This system 
provides a hard upper limit on the stellar mass inside its Einstein radius, no matter the 
IMF model. If we assume that the SSP models are correct and that this galaxy has a high 
[Na/Fe] abundance, we infer an IMF with a power-law slope x = 3-3.5 (where the IMF 
follows dn/dm = m~ x , and the Salpeter slope is x = 2.35). To assess whether these steep 
IMF slopes are consistent with the upper limit on the total mass, we calculate the total 
luminosity and the SSP stellar M/L ratio in stars for each assumed IMF to infer the stellar 
mass fraction inside the Einstein radius (-R E m — 4.55 ± 0.23 kpc; Koopmans et al. 2006). 
Changes in the IMF of stars with M < 0.3 M changes the total luminosity of the lens 
galaxy by at most ~ 10%. Conversely, stars with masses of 0.1-0.3 M & contribute > 60% of 
the stellar mass for bottom- heavy IMFs with slopes steeper than Salpeter (see, e.g., Fig. 2 of 
CvD12). To determine the stellar M/L ratio, we use the isochrones at solar [Fe/H] and [ct/H] 
for a 13.5 Gyr population from the Dartmouth Stellar Evolution Program (DSEP), a state- 
of-the-art stellar evolution code (Chaboyer, Green, & Liebert 1999; Chaboyer et al. 2001). 
We compare three different IMFs: Salpeter (x = 2.35), a bottom-heavy IMF (x = 3.0) and 
a very bottom-heavy IMF (x = 3.5). CvD12 use the same isochrones in their SSP for the 
bulk of the main sequence and red giant branch, except at M < 0.2 M , where they use the 
Baraffe et al. (1998) isochrones. For each IMF we compute the quantity 

= M*/M Ein = (L Ein /M Ein ) x (M*/L) DSEP , 

where M Ein is a robust measurement of the total mass enclosed within the physical Einstein 
radius [M Ein = (39.6 ± 0.8) x 1O 1O M ], L Ein is the luminosity enclosed within the Einstein 
radius, evaluated using B-spline luminosity models, as a fraction of de Vaucouleurs total 
model luminosity [L Ein = (4.49±0.2) x 10 10 L Q , from Bolton et al. 2008], and (M*/L) DSEP is 
the mass-to-light ratio from the DSEP isochrone using the appropriate IMF. The stellar M/L 
ratio includes the contribution from stellar remnants and gas ejected from stars at the end 
of their life-cycles. We list the results of this calculation in Table [2] For a Salpeter IMF, the 
stellar mass fraction of SDSSJ0912+0029 in the restframe V-b&nd is / Ei n,Sai P = °- 59 ± °- 15 > 
in agreement with previous results (0.60 ± 0.09, Auger et al. 2009). The mass-to-light ratio 
calculated from the DSEP isochrone for a Salpeter IMF is M/Ly = 7.2 ± 2 (M/L) in the 
V band and M/Lg = 10.2 ± 3 (M/L) in the B band. The latter value is consistent with 
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the upper limit of M/Lb < 9.08 (M/L) derived from dynamical models of Barnabe at 
al. (2009) under the maximum bulge hypothesis. An IMF slope of x = 3.5 yields M/Ly = 
29±9 (M/L) , and M/L B = 43±13 (M/L) corresponding to /| inj3 . 5 = 2.4±0.8, inconsistent 
with the total lensing mass within the Einstein radius at the > 95% confidence level. An IMF 
of x = 3.0 in £>-band is also excluded at the > 90% level, as this corresponds to a fraction 
/Ein3 o b = 1-6 ±0.5. For both of the bottom heavy IMFs in 5-band and for the x = 3.5 IMF 
in V^-band, we obtain a stellar mass fraction within the Einstein radius in excess of unity, 
thereby violating the lensing constraint on the total mass of the system at the > 90% CL. 
The x = 3 model is only marginally consistent in \^-band, but /gin 3.0 v = 1-4 ± 0.4 implies 
that there is no dark matter within the Einstein radius. 



4.2. Systematic uncertainties 



The uncertainty on the value of f^ in has a number of contributions. The uncertainties 
in the mass and luminosity determinations from lensing are much smaller than differences 
in the values of M/L arising from the use of different stellar population evolution models. 
The emerging picture is that, for a fixed IMF, it is difficult to constrain M/L estimates to 
much higher accuracy than 0.1 dex (Gallazzi et al. 2008; Marchesini et al. 2009; Longhetti 
& Saracco 2009, Conroy et al. 2009, 2010). We examine mass-to-light ratios predicted for 
different IMF from different stellar population models in the rest-frame V- and -B-bands and 
compare predictions from Worthey (1994), Bruzual & Chariot (2003), Maraston (2005), and 
Vazdekis et al. (2010) for single stellar populations with ages 11.2-14.1 Gyr, solar (Z = 0.02) 
or super-solar metallicity [Z = 0.05). For each SSP and each IMF, we calculate an average 
value and a standard deviation that we associate with the inferred values of M/L. Changing 
the [Fe/H] abundance from to 0.22 yields a ~ 9% uncertainty on M/L, while changing the 
age of the stellar population changes Mj L by ~ 20% at fixed IMF. The latter is the dominant 
contribution to the final uncertainties. We propagated these errors into the stellar mass 
fraction. The errors on the stellar mass fractions in Table [2] include both the random error 
contribution and the systematic uncertainties due to the use of different set of isochrones, 
bands, and stellar population age and metallicity uncertainties. 



5. Conclusions 

In this letter we have studied the Na I A8190 and TiO A6230 features - both indicators 
of low-mass (< 0.3 M@) stars in massive ETGs - as a function of each other, of age and 
metallicity indicators (Mgb, Fe, H/3), of NaD, and of stellar velocity dispersion. We find the 
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following: (1) The observed Nal-NaD trend depends strongly on stellar velocity dispersion 
of ETGs and only match current state-of-the-art SSP models for ETGs with a < 250 km s^ 1 . 
The most extreme Nal index strength in our sample is found in a gravitational lens system, 
which should have an IMF slope x > 3 based on the best current SSP models. The total 
enclosed mass of this system, however, excludes slopes steeper than x = 3.0 at the > 90% CL 
or slopes steeper than x = 3.5 at the > 95% CL. We conclude that the NaD feature is still 
affected by as-of-yet not understood processes in the more massive ETGs (<r > 250 km s -1 ). 
A full spectral comparison, in combination with lensing and dynamical constraints, is planned 
to further strengthen these results and assess whether Nal and NaD (in some instances) are 
contaminated. (2) We find that the TiO feature at A ~ 6230A (Ti02) is a particularly 
promising feature to decouple the IMF from age, metallicity, and abundance pattern of the 
stellar population, especially when combined with metallicity-dependent indices. We find 
that this feature correlates well with Nal, if the two most extreme discussed in the 

text, are excluded. This correlation can be a crucial piece of evidence against interstellar 
contamination of the Na I A8190 sodium absorption lines, although this does not solve the 
problem of NaD absorption. If strong Nal features are indeed not due to ISM contamination, 
very massive ETGs have higher [Na/Fe] abundances (> 0.3 dex) and slightly bottom-heavy 
IMFs, correlated with their stellar velocity dispersions. We also find a clear trend of an 
increasing IMF slope between a = 200 to 335 km s -1 from Salpeter (x = 2.35) to x ~ 3.0, in 
agreement with the XLENS system, which excludes steeper IMFs at the high-mass end. 

Our results are the first SSP-based indications of a steepening of the low-mass end of the 
IMF with increasing galaxy mass within the class of LRG/ETGs. Our results (i) support a 



similar trend first found by Treu et al. (2010), (ii) extend the evidence based on SSP models 
that the IMF steepens from spiral to early-type galaxies (vDCIO), (iii) suggest that Nal and 
NaD (in some instances) could be contaminated by interstellar absorption, and (iv) support 



a similar trend found by Cappellari et al. (2012) based on stellar kinematics. The upper 



limit of x < 3, based on one of the most massive ETGs in our sample, a gravitational lens, 



also supports our previous similar finding (Spiniello et al. 2011). 
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